Topology of the polarization field in ferroelectric nanowires from first principles 
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The behaviour of the cross-sectional polarization field is explored for thin nanowires of barium 
titanate from first-principles calculations. Topological defects of different winding numbers have 
been obtained, beyond the known textures in ferroelectric nanostructures. They result from the 
inward accommodation of the polarization patterns imposed at the surface of the wire by surface 
and edge effects. Close to a topological defect the polarization field orients out of the basal plane in 
some cases, maintaining a close to constant magnitude, whereas it virtually vanishes in other cases. 



The drive towards ever smaller ferroelectric devices [1^ 
has resulted in novel geometries such as nanowires [2j |3] , 
nanotubes [4' "5", "F, and nanodots [7, 8 . It was ini- 
tially assumed that as the size of such nano-shapes con- 
tinued to decrease, the ferroelectric polarization would 
vanish: ferroelectricity is a collective phenomenon after 
all. However, seminal works by Fu, Bellaiche and Nau- 
mov [9l[T0] suggested that, before vanishing, the polariza- 
tion may first undergo a size-induced transition by which 
the dipoles form a ferroelectric vortex, a polar configu- 
ration that does not exist in bulk and persists down to 
very small diameters. This insight has triggered much 
work, both theoretical [IH [12 [13 HHTS, 16^ and ex- 
perimental [iTl HH Hn] towards establishing true critical 
diameters for ferroelectricity, understanding the charac- 
teristics of ferroelectric vortices, and establishing their 
existence. 

In the present work we have used ab-initio methods to 
calculate the polar configuration in ultra-thin nanowires 
of the archetypal ferroelectric BaTiOa (BTO). Our re- 
sults show that the topological landscape of the polar- 
ization field P(r) in nanowires is much more complex 
than hitherto assumed, with polar configurations other 
than vortices appearing. The new topological textures 
include saddle points and quadrupoles as well as more 
complex configurations, all of them well described within 
the classic framework of topological defects of different 
winding numbers [20 . Such patterns are induced by the 
surfaces and can be understood in terms of simple sur- 
face and edge effects. Of equal importance is the discov- 
ery that the axial component of the polarization can ei- 
ther decrease or increase towards the surface of the wire 
depending on its termination, thus resolving the long- 
standing controversy on the sign of the so-called extrap- 
olation length [211 [111 [IS] . More importantly, the critical 
size for the persistence of polarization in nanowires is 
found to be determined by their surfaces. 

The two-dimensional (2D) topology of of the polariza- 
tion field in thin nanostructures can be mapped onto the 
topology of {Px^Py) for wires if the field remains con- 
stant along the axial direction of the wire. This paper 
focuses on topological defects in such 2D field [20 (these 



so-called "defects" do not refer to structural defects of 
the wires). The winding number is defined as n = 0/27r, 
being (j) the total angle the 2D field vector rotates when 
going around a closed circuit. For a continuous field n 
must be an integer, which if different from zero implies 
the existence of at least one defect within the circuit. 
Any smooth deformation of a vector field will conserve n, 
while changing the winding number implies a dramatic, 
extended and highly energetic rearrangement [33^. Such 
An 7^ rearrangements have been used to define "phase 
transformations" [24] for nano-objects displaying vortices 
(n = +1). Here we find realizations ofn = 0, 1, —1, and 
—3. We see textures that were described before (e.g. ra- 
dial [25]) together with completely novel ones. The differ- 
ent textures are induced (and thus can be manipulated) 
by the different surface terminations of the wires. 

Our calculations were performed within DFT and the 
generalized gradient approximation (GGA) [26^. We used 
the Siesta method [27\ [^, based on finite-range nu- 
merical atomic orbitals, using a double-^ polarized ba- 
sis set [m [29] . Norm-conserving pseudopotentials were 
used, including into the valence the semicore shells 35 and 
3p for Ti and 5s and 5p for Ba. The performance of the 
method was tested for bulk BaTiOs. The obtained lattice 
parameter for the rhombohedral phase is 0.12 % smaller 
than the experimental value of 4.0036 A at 15 K [30], and 
0.055 % smaller than the published results with the same 
GGA and plane waves as basis set [26]. The rhombohe- 
dral angle was found to be 89.84°, exactly as obtained 
in the same experiments, and the off-centering of the Ti 
atom came out as 0.184 A, as compared with the experi- 
mental 0.185 A. Periodic boundary conditions were used 
with a simulation cell describing a square array of infinite 
nanowires, such that periodic replicas are separated by 
24 A. The width of the wires is shown in Fig. 1. The unit 
cell along the wire axis comprises a single bilayer in all 
cases, and therefore all polarization textures are homo- 
geneous along z by construction. Integrals in reciprocal 
space used a ^-mesh of 10 A cutoff [31 , while the inte- 
grals in real space used a r-mesh of 350 Ry cutoff [28 . 
The atomic positions were relaxed until all atomic force 
components were smaller than 10 meV/A. From the re- 
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FIG. 1: (Color online) Sections of the different BaTiOs 
nanowires considered. The largest balls (cyan) are for Ba 
ions, the intermediate ones (magenta) represent Ti ions, and 
the smallest ones (red) are for O. (a) is Ti02 terminated, (b) 
is BaO terminated, except the edges, where lines of Ba atoms 
have been removed, (c) is stoichiometric, with two BaO sur- 
faces and two Ti02 ones, (d) is BaO terminated. 



laxed structures, the Ti displacement from the centre of 
its coordination octahedron was taken to define P(rTi), 
as a measure of the local dipole. For Ti atoms on surfaces 
or edges the off-center displacement was determined by 
the remaining coordinating O atoms. It should be em- 
phasized, however, that the focus of the paper is on the 
topological aspects of the polarization, which are not af- 
fected by quantitative details related to this definition. 

Four different kinds of nanowires have been considered, 
corresponding to different surface terminations, as shown 
in Fig. 1, all showing (100) and (010) surfaces, or slight 
alterations thereof. Two of them display the same termi- 
nation on all sides, namely, all Ti02-terminated (a), or 
all BaO-terminated (d). Wire- type (c) is stoichiometric, 
with both terminations, and corresponding to the wires 
studied in ref. [11 . Finally (b) is like (d) except for the 
fact that the Ba on the edges has been removed. Each 
wire type has been considered in several sizes, mainly 
4x4, and 5x5, counting the Ti atoms. 

Fig. 2 shows the polarization textures of the wire types 
of Fig. 1, alongside qualitative sketches of the correspond- 
ing topological defects for a 2D vector field. The discrete 
set of off-centering vectors is assimilated to a lattice dis- 
cretization of a 2D P field. Different wire widths cor- 
respond to different discretizations of the corresponding 
general pattern. A first striking result is that even for the 
small sizes considered (and the coarse discretization they 
imply) the field textures appear very clearly. All textures 



preserve inversion symmetry at the center of the wire 
if it was not previously broken by the surface termina- 
tion. The first two cases in Fig. 2, displaying inward and 
outward radial patterns, correspond to the same kind of 
topological defect, of winding number 1. They are thus 
homotopical with each other, but also with the vortex 
structures seen in refs. [TOl [13]. The continuous trans- 
formation taking a radial texture onto a vortex one is 
beautifully illustrated in Mermin's review [20]. They of 
course differ in other physical aspects, such as the non- 
zero toroidal moment of the vortex patterns. 

Figs. 2(c) and (d) correspond to different winding num- 
bers as yet not observed or proposed in this field, namely 
— 1 and —3. In Fig. 2(c) the defect is displaced from 
the center (it is the only non-centrosymmetric wire). 
Fig. 2(d) shows the n = — 3 case of the BaO terminated 
wire. The high winding number imposes a more rapidly 
varying field around the defect, which implies higher en- 
ergy. The system responds by a more pronounced depres- 
sion of the magnitude of the 2D field around the defect. 
The winding of the field around the outer circuit in the 
wire section clearly corresponds to n = —3. 

Fig. 2(e) is for the same BaO termination as Fig. 2(d) 
but smaller thickness. Although a different discretiza- 
tion, the winding around the outer circuit clearly shows 
the same n = — 3. In this case the inner displacements, 
although small, all point mainly outwards, inconsistently 
with the the n = — 3 texture: This system has found 
more stability in decomposing the global n = — 3 global 
defect into a set of defects whose winding numbers add 
up to —3. In general [20 , when following the winding of 
the field around a circuit the winding number obtained 
is the sum of winding numbers of all the point defects 
enclosed in the circuit. In this case the texture can be 
described by a central defect with n = +1, the field radi- 
ating outward as in Fig. 2(g), plus four other defects of 
n = — 1, as illustrated in the last sketch. 

The described textures are clearly induced by surface 
and edge effects. Indeed, the phenomenlogy described 
can be understood as originated by three simple ten- 
dencies: {i) In Ti terminated surfaces and edges, the Ti 
atoms tend to off-center inwards, (n) Ba surface termi- 
nation induces outward Ti off-centering patterns, while 
{Hi) Ba edge termination induces an inward tendency. 
The competition of (ii) and {Hi) in the Ba terminated 
wires gives rise to the rich n = — 3 pattern. By remov- 
ing the Ba edges [Fig. 2(b)], the competition disappears 
and the fully outward pattern is established. Similarly, 
when removing the Ba edge in the stoichiometric case 
the polarization field becomes topologically homogeneous 
(n = 0), with the field lines crossing the wire diagonally 
pointing left and down (not shown). Special chemical 
tendencies at edges and surfaces have been previously 
reported [161 HZj- Shimada et al. [16] talks specifically 
about strong "edge bonds" for PbTiOs nanowires. 

The surfaces not only determine the (P^^, Py) topology. 



3 



(a) 



\ \ A A 



(b) 



(c) 



(d) 



(e) 



\ 



i (f) 



(g) 





1 


1 


/ 


/ _ 


N 




/ 


/ 








/ 


/ 




\ 




\ 


1 




/ 




\ 


\ 





\ 




f 






/ 








\ 












\ 








/ 


/ 


V 




y 


\ 



\ 


/ 


\ 


/ 








V 










/ 


\ 


/ 


\ 




FIG. 2: Left column (a)-(d): Ti off-center displacements for 
the four respective nanowires in Fig. 1. (e) is analogous to 
(d) but for a thinner wire. The right column sketches the cor- 
responding 2D field lines around topological point defects of 
winding numbers n = 1, 1, — 1, and —3, respectively, (j) show- 
ing the n = — 3 decomposition into a central n = +1 and four 
n——\ defects. The respective largest arrows corresponds to 
(a) 21 pm, (b) 35 pm, (c) 23 pm, (d) 6.8 pm, and (e) 6.6 pm. 



but also the behaviour (Fig. 3). Axial switchable po- 
larization has been observed experimentally p^, and has 
been analyzed within Landau- Ginzburg theory p!2| [T4]. 
and DFT [HI [16]. We find that for TiOs termination 
the 4x4 wires display = 0, while it is finite for 5 x 5 
(in the figure) or larger. When non-zero, the axial polar- 



ization stems mainly from the surface cells, the central 
cells displaying much smaller displacements. This is sur- 
prising: an essentially constant-magnitude field is com- 
monly assumed for these and other ferroelectric systems 
(some models have this assumption built in, as a constant 
dipole moment of arbitrary orientation). Here, however, 
the field responds by vanishing smoothly into the defect. 
The BaO terminated wires respond in the more expected 
way: discontinuity at the defect is avoided without sub- 
stantial diminishing of the polarization field, by having 
maximal axial polarization at the wire center. This is the 
case for both with and without Ba edges, although the 
latter shows a much smaller axial component at the cen- 
ter. There is also a 3 x 3 critical thickness for (largest 
system with P^ — 0) in the Ba terminated case. In the 
stoichiometric wires the situation is different. The figure 
shows a 5 X 5 case (the critical thickness for this type 
is 2 X 2, in agreement with ref. [11 ), which again shows 
large P^ for the surface and edge Ti displacements, for 
which the in plane components were also large. 

In addition to their instrinsic fundamental interest, the 
polarization-field textures and the possibility of manip- 
ulating them, offer interesting exploitation possibilities. 
Net polarization has been measured and switched per- 
pendicularly to the wire axis at room temperature with 
a scanning probe microscope [U [3l [17]. Such a setup 
should thus respond as a three-state system, with +1,0 
and -1 states, the +1 and -1 corresponding to the conven- 
tional ferroelectric response, and the zero to the stable 
topologically defective state (with \P\ much smaller than 
for the other states, strictly zero if the unpolarized wire is 
centrosymmetric). The field rearrangement needed for a 
non-homotopical change of texture would mean an energy 
barrier when changing the net polarization, giving a sce- 
nario analogous to a first order phase transition above the 
transition temperature: the global minimum for P ~ 
and two local minima for + and -. Note that some of 
the n 7^ textures exhibit net quadrupoles. These con- 
figurations can be switched (destabilized) by application 
of an electric field gradient, which could in principle be 
used for a quadrupole memory device. 

Phase transformations were already proposed for vor- 
tex structures in nanodots [10 . In wires, as extended 
systems, the expression would be more appropriate ex- 
cept for the fact that wires with sufficiently short-ranged 
interactions do not sustain spontaneous symmetry break- 
ing at any finite temperature [32^. Even assuming that 
the long range dipole-dipole interactions and the also 
long-ranged strain mediated interactions do not invali- 
date the previous argument, the discussions in this pa- 
per can be understood in terms of the observed wires 
being of lengths smaller that the dipole-dipole correla- 
tion length: any perturbation pins the order parameter 
(the polarization across wires is experimentally observed 
after all). As a final remark, the surface and edge control 
not only affects the polarization across the wires, but can 
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FIG. 3: z component of the Ti off-centering for the wires in 
Fig. 1. The largest circle represents a value of 15 pm. 



also determine the conductance along the wire, since the 
deviations from (BaO)^E(Ti02)y compositions give rise to 
doped wires, opening still further possibilities for these 
systems, to be explored elsewhere. Evidently, different 
surface chemistries than the ones contemplated here can 
be achieved depending on how the wires are produced. 
The points remain, however, that interesting topologies 
are to be expected in thin ferroelectric wires, and that 
they are induced by the surfaces and edges. 

In summary, based on first-principles calculations, we 
predict ferroelectric patterns in BTO nanowires for differ- 
ent surface terminations, finding new textures that corre- 
spond to the topological winding numbers n = 0, +1, — 1 
and —3 of point defects in 2D vector fields. A tendency of 
the n = —3 defect to decompose into one +1 and four —1 
defects is observed, as well as the different mechanisms 
used by the field to avoid discontinutities. 
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arise for vector fields of constant magnitude, giving rise to 
a discontinuity in the field derivatives at the defect. Ferro- 
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are of variable magnitude and vanish at the defect. 



